Field Evolution of Coexisting Superconducting and Magnetic Orders in CeCoIn 5 
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We present nuclear magnetic resonance (NMR) measurements on the three distinct In sites of 
CeCoIns with magnetic field applied in the [100] direction. We identify the microscopic nature of 
the long range magnetic order (LRO) stabilized at low temperatures in fields above 10.2 T while still 
in the superconducting (SC) state. We infer that the ordered moment is oriented along the c-axis 
and map its field evolution. The study of the field dependence of the NMR shift for the different 
In sites indicates that the LRO likely coexists with a modulated SC phase, possibly that predicted 
by Fulde, Ferrell, Larkin, and Ovchinnikov. Furthermore, we discern a field region dominated by 
strong spin fluctuations where static LRO is absent and propose a revised phase diagram. 



The detrimental effect of an applied magnetic field on 
a superconductor has nourished the intuitive presump- 
tion of a de facto competition between magnetic and 
superconducting (SC) orders. However, it is now well- 
established, both theoretically and experimentally, that 
not only can these orders coexist, but in some cases, 
they may even be essential for each other's stability [![. 
Manifestations of coexistence span a rather wide range of 
materials including several cuprates, ferropnictides, and 
heavy fermion systems of which CeCoIn 5 is one of the 
most intriguing examples Q ■ In the SC state of this com- 
pound application of a magnetic field (Ho) induces a long 
range magnetic order (LRO) , restricted to a narrow low- 
temperature (T) region of the phase diagram below H c i 
[H, 01 ■ What is more, this particular region of the phase 
diagram was initially identified as the first realization of 
the long-sought Fulde, Ferrell, Larkin, and Ovchinnikov 
(FFLO) state, a superconducting state with a non-zero 
pair momentum and a spatially modulated order param- 
eter 0, @|. However, important questions regarding the 
true nature of the low-T high-i/ SC phase, the details 
of the magnetic order and its field dependence, and the 
potential driving mechanisms of their coexistence remain 
unanswered 0]. Thus, CeCoIn 5 provides a strikingly rich 
ground to study the complex interplay between exotic SC 
and magnetism. Experimentally, nuclear magnetic reso- 
nance (NMR), as a microscopic probe sensitive to both 
magnetic and SC degrees of freedom, provides a powerful 
tool for the investigation of these puzzles. 

In this letter, detailed low temperature NMR mea- 
surements on the three distinct In sites in CeCoIns for 
H 1 1 [100] are presented. We establish that at T w 70 mK 
a phase with static magnetic LRO is stabilized for fields 
above w 10.2 T in the SC state. We deduce that the LRO 
is an incommensurate spin density wave (IC-SDW) with 
moments oriented along the c-axis, independent of the 
in-plane Ho orientation. Further, the detailed field evo- 
lution of the moment is mapped. The study of the field 
dependence of the NMR shift implies that this IC-SDW 



coexists with a novel SC state, characterized by an en- 
hanced spin susceptibility Q. Finally, we identify a new 
region in the Hq-T phase diagram, lying in between the 
low field SC (lfSC) and the IC-SDW states. This could 
be an FFLO phase without magnetic LRO, in agreement 
with recent theoretical predictions [9(. 

High quality single crystals of CeCoIn 5 , grown by a 
flux method, were placed in NMR radio frequency (RF) 
coils inside the mixing chamber of a dilution refrigera- 
tor so that Ho || [100]. The RF coil was used to deter- 
mine the precise value of Hq by performing 63 Cu NMR 
on its copper nuclei. The spectra were obtained, at each 
given value of Hq, from the sum of spin-echo Fourier 
transforms recorded at constant frequency intervals. Ex- 
tremely weak RF excitation power [1 01 ] was used to dis- 
cern the NMR signal from In sites sensitive to magnetism. 

For Ho || [100], there are three inequivalent In sites. 
The axially symmetric In(l) is located in the center of 
the tetragonal Ce planes, while In(2 ac ) and In(2f, c ) sites 
correspond to In atoms located on the lateral faces (paral- 
lel and perpendicular to the applied field, respectively) of 
the unit cell [ll|. In Fig. [T]the H evolution of the In(2 ac ) 
and In(2;, c ) spectra at T ~70 mK is plotted. Lowering 
the field below ~ 11.7 T establishes magnetic LRO. The 
LRO is evident in the fact that the In(2 ac ) line broad- 
ens into a spectrum with two extrema/peaks with finite 
signal weight in between them. Such spectra are charac- 
teristic of IC LRO along one spatial dimension [12j . Fur- 
thermore, at Ho = 11.67 T, both the broad In(2 ac ) and 
the sharp normal state spectra are observed. This reflects 
the coexistence of normal and IC phases in the vicinity 
of the phase transition confirming its first order charac- 
ter. From independent NMR measurements [3] and the 
tuning resonance of the tank circuit, we establish that 
this transition to the IC LRO state coincides precisely 
with the transition from the normal to the SC state. For 
9.2 T < H < 10.2 T, the spectra of all In sites consist of 
a single peak, i.e. no signature of the IC state is observed. 
However, these spectra remain significantly broader than 
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FIG. 1: (Color online) NMR spectra of (a) In(l), In(2 ac ) 
and (b) In(2 bc ) at T = 70 mK for various H || a. The 
frequency scale is denned by subtracting u>q, the zero NMR 
shift frequency. N denotes the normal phase, IC the LRO 
phase, eSC the state with strong fluctuations, and lfSC 
the Abrikosov SC state. Solid lines in (a) are simulated 
spectra for the IC-SDW order described in the text. 

the ones for Ho < 9.2 T, where the linewidth of all sites 
can be adequately described by the spatial distribution 
of magnetic fields resulting from the vortex lattice [l3| ■ 

We now proceed to the analysis of the nature of the IC 
LRO phase to elucidate the field evolution of its magnetic 
moment and its ordering wave vector (Q). The NMR 
spectrum reflects V(Hf nt ), the probability distribution of 
the internal magnetic field at the nuclear site projected 
along the Ho direction, denoted by z. In the LRO phase, 
the component of the internal hyperfine field parallel to 
H , at an In site, is given by Hf nt = z ■ ■ Mi, 

where A,; is the symmetric 3x3 hyperfine coupling ten- 



sor with the i nearest-neighbor Ce atom and fii is its 
magnetic moment 
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The exact form of the hyper- 
fine tensor is derived in Ref. [l5|. Our data and the 
fact that certain off-diagonal elements of A (A a b, A ac , 
and Ab c for In(l), In(2 ac ), and In(2& c ), respectively) are 
non-zero allows us to place stringent constraints on the 



possible nature of the magnetic LRO. Specifically, the 
form of Q and fii should be such that the In(2 ac ) spec- 
trum (Fig. [TJi) broadens into a double peak structure, 
while no such broadening of the In(l) [l3j and In(2(, c ) 
(Fig. [T]d) lines is induced. In principle, such spectra can 
be effectively described by an appropriate sinusoidal vari- 
ation of -ffJL, that is, of the moment itself. For an IC- 
SDW state, this spatial variation of the moment can be 
written as (J,i(r) — fio cos(Q • r + fa) where Q is the IC 
wave vector, r defines the lattice coordinates of the i th 
Ce moment, and fa is an arbitrary phase. The separa- 
tion between the two extrema of the In(2 ac ) spectrum 
is Af = "f(H\\^ ax — H^" 1 ), where 7 is the nuclear gy- 
romagnetic ratio. Thus, Af depends on a product of 
hyperfine tensor components, the moment amplitude /io, 
and some trigonometric function dictated by Q. 

A priori, the magnetic structure with the ordered mo- 
ments /jj||Ho||a and Q = Qaf + 8, where Qaf = 
(0.5,0.5,0.5) and incommensuration 8 = (5,0,0), can 
satisfy the requirements imposed by our data, as was 
suggested in [3j. In this case, Af for the In(2 ac ) peaks 
is equal to 4jA aa IM) sin(^7r), with A aa w AfJ^^/2 w 
— 0.6T//1B (TH. flB]]. However, reproducing the observed 
spectra requires cither a moment fiQ that is three times 
larger than that found by neutrons for Ho|| [110] |H or an 
excessive value of 8 ~ 0.21. It is thus unlikely that such 
a structure with fj,i\\a is stabilized in the IC phase. 

A better candidate for describing our data is an IC- 
SDW with ordered moments perpendicular to the plane, 
i.e. Hi\\c L H as in and 8 in the plane. We cannot 
uniquely determine the direction of 8, since the relevant 
V(Hf nt ) is nearly the same for <$||a, 8\\b, or <$||[110]. How- 
ever, the case of 8\\b _L Ho appears to be the most plausi- 
ble one, since it would agree with both the experimental 
findings for Ho || [110] [4j and the theoretical prediction 
of a 8 _L Hp for coexisting IC-SDW and FFLO order pa- 
rameters [9|. In this case, Af is proportional to 4rfA ac ^, 
where A ac is the off-diagonal hyperfine tensor compo- 
nent, whose value is not known from other independent 
measurements. Assuming that at 11 T un = 0.15/^b and 
\S\ = 0.085 (the same magnitude as in [J]) we find that 
a reasonable value of A ac « 0.38 T//ig is required to 
fully account for our data. The calculated distributions 
V(Hf nt ), convolved with the underlying vortex lattice 
lineshape are depicted as the solid lines in Fig. QJi. 
Based on this result, the magnetic LRO for H | \a is most 
likely an IC-SDW with the ordered moment perpendicu- 
lar to theplane. Hence, in conjunction with the neutron 
findings [4| , we infer that the direction of the ordered mo- 
ments is independent of the in-plane orientation of Hq. 

Having inferred that the magnetic moment is ori- 
ented perpendicular to the plane, we proceed to extract 
the field evolution of its magnitude /j,q from Af of the 
In(2 ac ) depicted in Fig. [Ik,. To do so, we fit the data to 
V(Hf nt ) with /ip as a fitting parameter and assuming that 
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FIG. 2: (Color online) Magnetic moment as a function of 
Ho deduced from A/ at T » 70 mK. The double arrow 
denotes the region of Ho within which LRO is absent but 
spectra remain much broader than in the lfSC state. 

\S\ = 0.085 is Ho independent, as exp erimentally found 
Q and theoretically predicted @, Il6| . and that the hy- 
pcrfinc component A ac « 0.38T//is. The deduced field 
evolution of fiQ is shown in Fig. [2] The moment is zero 
outside the IC phase below Hq 10.2 T. It increases by 
a factor of ~ 2.5 as the field changes from « 10.3 T to 
11.67 T, where it reaches its maximum value of Fts 

In order to better understand the perplexing relation- 
ship between the coexisting SC and IC-SDW orders, we 
next discuss the field dependence of the NMR shift of 
In(2 ac ) and compare it to that of In(l) [l3|. This al- 
lows us to distinguish between the contribution from the 
low energy local density of states (LDOS), characteristic 
of the SC state, which affects both sites and that from 
the localized magnetic moments, affecting predominantly 
In(2 ac ). Our previous analysis of the In(l) shift has re- 
vealed that a novel SC state, possibly an FFLO state, 
is stabilized for Hq = H* > 10 T via a 2 nd order phase 
transition [13| . We were, however, unable to exclude the 
existence of a magnetic order in this novel SC phase. In 
Fig. [31 the field dependence of the shift, determined by 
diagonalizing the full nuclear spin Hamiltonian, is plot- 
ted. In the lfSC phase, the DOS increases with Hq due 
to excess Zeeman and Doppler-shifted nodal quasiparti- 
cles [HI], leading to an increasing spin susceptibility 
The In(l) and In(2 ac ) shifts (K cx A\ s ) exhibit reversed 
Hq evolution. That is, they scale with their respective 
hyperfine coupling constants, which are nearly equal but 
have opposite signs with that for In(2 ac ) being neg- 
ative. The same is also true for the relative shift change 
between the normal and lfSC state for both sites. 

For H* < H < H c2 , where the IC-SDW LRO is es- 
tablished, we plot the In(2 ac ) shift of its first moment 
and of the two extrema of its broadened lincshapc. The 
difference of the shift of the two extrema strongly varies 
with the field reflecting the increase of A/ with increas- 
ing Hq. Further, the In(2 ac ) shift associated with its first 
moment has essentially the same field dependence as that 
reported for In(l), once the difference in their hyperfine 
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FIG. 3: (Color online) NMR shift (spin, K 3 , and constant 
orbital part, K orb [ll[) of In(l), In(2 ac ) and ln(2 bc ) at 
T ~ 70 mK as a function of Ho||fi. The solid lines indi- 
cate the values of H* and H C 2- 

couplings is taken into account. Thus, independently of 
the local magnetism, both In sites sense an additional 
common local field, due to a spin polarization along Ho . 

There are two likely origins for this local field for 
Hq > H* . One is a canting of the transverse staggered 
magnetization of the localized moments along Ho. To 
test this hypothesis one needs a quantitative description 
of canting in an IC phase, which is missing. Thus, we 
calculate the effect of several crude models of canting 
on the NMR observablcs. We find that our data are ef- 
fectively reproduced by considering an appropriate fixed 
value of /Xf (~ /io sm (5 )) on all lattice sites, regardless 
of |/n,(r)|. Nevertheless, it is not clear how this type of 
canted moment could be induced in an IC-SDW phase. 

Alternatively, the additional common field could orig- 
inate from the enhanced LDOS of the spin polarized 
quasiparticles in the nodal planes of the FFLO state [|| . 
In that case, since in an FFLO state LDOS varies on 
a length scale that significantly exceeds the spacing be- 
tween two In sites, its contribution to the shift of the two 
distinct In sites should be equivalent as observed. This 
spatial modulation of the LDOS can also contribute to 
asymmetric spectral broadening, which for different In 
sites should scale as their respective hyperfine couplings. 
For In(l), this LDOS modulation gives rise to a tail on 
the high frequency side of the line [l(3, [l3j], while for 
In(2 ac ) the tail should be on the low frequency side, due 
to the sign difference of their hyperfine coupling. This is 
indeed observed for In(2 ac ) spectra, where only the low 
frequency side is essentially broadened beyond IC-SDW 
lincshape as shown in Fig. 

Next, we consider the intriguing low field limit of the 
IC LRO phase. The low-T shift data provide evidence of 
a continuous phase transition at « 10.2 T, as shown in 
Fig. [3J However, significant broadening of ln(2 ac ) and 
In(2b c ) spectra, as compared to the ones in the lfSC 
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FIG. 4: (Color online) A sketch of an alternative phase 
diagram of CeCoIns (data points adopted from Ref. P3). 
The solid and dashed lines indicate 1 st and 2 nd order phase 
transitions, respectively. 

phase, onsets for Hq > 9.2 T, as evident in Fig. [T] Strik- 
ingly, 9.2 T is precisely the field at which the amplitude 
of the IC-SDW should vanish if one extrapolates from its 
field dependence in Fig. H For 9.2 T < H Q < 10.2 T, the 
In(2 ac ) signal is weak due to the rapid loss of spin coher- 
ence, caused by strong field fluctuations. Thus, although 
the static LRO is absent, strong antiferromagnetic (AF) 
fluctuations are still present. These fluctuations can be 
responsible for the apparent collapse, evident in the dis- 
appearance of the double-horned broad In(2 ac ) line, of 
the LRO in this field region. Besides, the observed addi- 
tional static broadening, as compared to the IfSC state, in 
the presence of these fluctuations, can only emerge from 
highly enhanced \s around defects [Tjj, such as nodal 
planes in an FFLO state. Thus, this field range possibly 
corresponds to an FFLO state in the presence of strong 
AF fluctuations, as predicted in [gj]. This is all more likely 
since the width of the low-T spectra of all three In sites 
in the same field range is comparable to that observed in 
the higher T (T > T*), high H (Hn > H*) SC phase 
outside the limits of the LRO state [H, [T3| ■ 

The 2 nd order phase transition previously identified 
by the In(l) shift can be, then, the transition to the IC 
LRO state. As recently shown, the IC magnetism can 
arise in the FFLO state in a <i-wave SC as a consequence 
of the formation of Andreev bound states near the zeros 
of the FFLO order parameter It is a large LDOS in 
the bound states that triggers the formation of the IC 
LRO . At finite T the LRO phase is stabilized only 
when a sufficient number of the nodal planes containing 
the bound states is induced by Hq jjj [l9| ■ Since these 
same states contribute to the shift, it is likely that the 
transition identified by it is indeed the one to the IC LRO 
phase. Alternatively, if the visible increase in the In(l) 
shift as compared to that in the IfSC phase were due to 
canting, an evident increase should onset only after LRO 
is established. 

Based on these observations and consistent with the 
discussion in the previous paragraphs, we postulate that 
a phase with strong AF fluctuations (referred to as 'ex- 
otic' SC (eSC)) develops for 9.2 T < H Q < 10.2 T at 



low-T, while the 2 nd order phase transition previously 
reported @,[H,[20[ at H* w 10.2 T marks the transition 
to the state with well established static IC-SDW order. 
Such a two step phase transition from IfSC to FFLO (as- 
suming it exists in eSC) and then to the IC state was 
theoretically predicted when spin fluctuations are consid- 
ered in Ref. [9j. In light of our NMR results, we propose 
the revised phase diagram sketched in Fig. [4] 

In conclusion, our comprehensive low-T NMR data 
provide a clear picture of the field evolution of the IC 
magnetism, and the magnitude of the magnetic moment, 
confined within the high-7?o low-T phase of CeCoLi5. In 
conjunction with neutron scattering results, we deduce 
that the direction of the ordered moment does not de- 
pend on the in-plane applied field orientation. Our anal- 
ysis of the field dependence of the NMR shift on different 
In sites indicates that the IC-SDW order likely coexists 
with an FFLO state. Finally, we identify a novel phase, in 
the field regime in between the IfSC and IC-SDW states, 
which is likely a true FFLO phase in the presence of 
strong AF fluctuations. 
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